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Abstract 

A simple microscopic model of a small superconducting loop interrupted by Josephson junction (flux 
qubit) allows to compute from the experimental data of Wal et.al [l| an important parameter - the density 
of Cooper pairs at zero temperature. This density is determined by the cut-off energy in the BCS model 
and agrees with the original BCS suggestion but is lower by two orders of magnitude than the value 
accepted in the modern literature. The immediate consequences of this result are: the validity of the 
strong coupling BCS model, a plausible picture of electrons recombination into Cooper pairs, and a much 
weaker condition for the appearance of high-temperature superconductivity. Another consequence is that 
the popular interpretation of Josephson qubits as macroscopic quantum systems is replaced by a picture 
of qubit states being superpositions of the ground state and the state containing only a single excited 
Cooper pair. 

PACS numbers: 74.20.Fg, 85.25. Cp 
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Fifty years after the discovery of superconductivity Bardeen, Schrieffer, and Cooper (BCS) 
presented a theory of this phenomenon which remained valid for the next fifty years |2j. The BCS 
theory contains a cut-off energy parameter determining the number of electrons participating in 
the pairing mechanism and initially suggested to be comparable to the superconducting gap. For 
many years this parameter remained unspecified [3] but in the modern literature has been replaced 
by the Debye energy Q wh|ch is about hundred tunes higher than the gap. The relatively new 
experiment with flux q ubit ,1, sheds a new light on that issue supporting the original hypothesis 
The detailed analysis of the microscopic Hamiltonians based on this hypothesis is presented in [5] 
for all three types of Josephson qubits (charge, phase, and flux qubits 6|]). Here, only a simplified 
model of flux qubit is discussed. 

A flux qubit (FQ) is a small superconducting ring interrupted by one or several Josephson 
junctions. The thickness of the ring is assumed to be not larger than the penetration depth for 
the magnetic field (typically ~ lOOnm), what allows to ignore spatial variations of the electric 
current density. The Cooper pair's circular states \k,fi), which produce current, are labeled by 
two quantum numbers. The first one k — 1, 2, .., K describes a Cooper pair in its center of motion 
reference frame and corresponds to a pair of time-reversed single electron states \k±). Under 
standard half filling condition there are K/2 Cooper pairs in the sample at zero temperature. 
The second quantum number /i = 0, ±1, ±2, ... accounts for the quantized circular motion of the 
Cooper pair along the ring. According to the Onsager hypothesis [7] a single Cooper pair in a 
state \k,/j.) generates a quantized magnetic flux /i$o where $o = h/2e. Consider a configuration 
of Cooper pairs described by the occupation numbers which take values 0, 1, as Cooper pairs 
behave like hard core bosons. The total magnetic flux is equal to F $ 0; where F = Ylk nt in kv an< ^ 
the magnetic energy of such configuration is given by 

£m = E L (^2fin kf ^ . (1) 

where 

E L = Kh 2 /Am£ 2 , (2) 

m is electron's mass and £ is the length of the loop. The formulas ((T|) and (|2]) are derived using 
the London's relation between superconducting current density and vector potential j = — ^-A, 
the formula for the energy of a current in a magnetic field £ — — J j • A d 3 x. = ^-A 2 and the 
Onsager's magnetic flux quantization condition $ = § A • c?x = ±A£ = -F$ , F = 0, ±1, ±2, ... 
Here, J\f is a density of superconducting electrons and A denotes |A|. Both quantities are assumed 
to be uniform in the sample what is true for thin enough loops (thickness < lOOnm). 



For loops of micrometer diameters, used to construct flux qubits, El is much higher than the 
superconducting gap what means that the states with different fluxes F $ are separated by large 
energy gaps in comparison with the other relevant energy scales for the system. To produce a qubit 
one has to switch on an external flux <& ext = fi ex t^o which modifies the magnetic energy yielding 

£m(/w) = E L (F — fl ext ) . (3) 

Taking fi ext ~ 1/2 one can reduce the magnetic energy difference 6~£m = Ej\ \ — 2^ ext \ between the 
unique ground state |0) with a flux F = and the excited levels with F — 1. The stable states of a 
superconducting ring correspond to collective states of Cooper pairs with various flux population 
numbers = Ylk n kv sucn th & t; Y^u^i^ = K/2. In particular the ground state |0) has the total 
flux equal to zero with concentrated around /i = 0. The stable states with F = 1 are obtained 
from the ground state by changing the populations N„ i— > N v — 1, N u+ i i-> N v+ i + 1 at the given 
v. Such states are denoted by |0,z/). The presence of Josephson junction(s) does not produce 
transitions between |0) and |0, v) because those states have a global and collective character while 
the Josephson junction should be described by a localized Hamiltonian perturbation Vj. However, 
the transitions exist between |0) and the states |0, u, k) obtained from |0, v) by replacing one ground 
pair with the indices [y + 1, k) by an excited pair. Excited Cooper pairs appear already in the 
original BCS paper [2] and possess a clear interpretation in terms of permutational symmetry within 
a strong-coupling BCS model The corresponding transition amplitudes £ uk = (0, v, k\Vj\0) 
determine the state 

\l) =Y,U\0^,k) , i,* = -=h=r= (4) 

u,k V l^v'k' IVfe'l 

which becomes separated from the others and together with the ground state |0) span the flux 
qubit Hilbert space. 

In the two level approximation the eigenstates of the qubit Hamiltonian are suitable superpo- 
sitions of |0) and |1) and the corresponding qubit frequency reads 



f = j-^[E L (l-2jl ext )Y + Ej. (5) 

Here, the Josephson energy Ej = ^2 vk \£, u k\ 2 and jx ext = fi ext + 5. The correction 5 accounts for the 
fact that the excited pairs have higher energy than the ground ones (see [5] for the details). 

The simple theory presented above can be now compared with the experiment of van der Wal 
et.al.{l|. The sample is an aluminum 5/im x 5/im loop made of 450nm wide and 80nm thick 
lines. The Josephson junction loop is kept at the milikelvin temperatures and excited by the 



microwave radiation with frequencies in the range of 1 — 10 GHz. The measured dependence of 
the resonant frequency on the control parameter jl ext in the range 0.495 < fi ext < 0.505 can be used 
to test the formula (jSJ). The off-set 5 is observed also, but attributed to persistent currents. The 
linear dependence on |1 — 2jl ext \ far enough from jx ext = 1/2 is clearly visible as well as the level 
repulsion for jx ext very close to 1/2 (see FIG.l). From those data the value E^/h ~ 1.5 x 10 3 GHz 




FIG. 1. Energy levels of the flux qubit (in GHz) as functions of the external magnetic flux (in 3>o units) 
given by (0). A fit to the data of van der Wal et.al. with E L /h = 1500GHz and Ej/h = OMGHz. 

is extracted. Putting t = 20/im into (j2J) one obtains K = 3.3 x 10 6 . The volume of the sample 
V = 0.72 x 10~ 12 cm 3 and hence, for the first time, one can compute directly from the experimental 
data the zero temperature density of Cooper pairs in a superconductor (Al) 

K OP [Al] = K/2V = 2.5 x 10 18 /cm 3 . (6) 

The obtained result can be discussed within the BCS theory 2|, |4j. The basic predictions of the 
BCS concern the relations between the superconducting gap at zero temperature A(0), the critical 
temperature T c , the coupling constant g describing the magnitude of electron-phonon interaction 



j^l, and the cut-off parameter hw c 

A(0) = 1.76 k B T c , k B T c = 1.13 - . . , . (7) 

2 sinn (nw c /g) 

The energy cut-off ftu; c limits the kinetic energy of the electrons participating in the BCS pairing 
mechanism to the interval [Ep — fad c , Ep + huj c ] (Ep - Fermi energy) and determines the number 
K of the corresponding electronic states by the expression 4] 

K = 2hw c N(0). (8) 

Here, N(0) is the density of electronic Bloch states (excluding electron's spin) at the Fermi surface 



given by 
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iV(0) = V 2^( 37r2 ^) 1/3 (9) 



where K e i is a density of electrons (K e ;[v4Z] = 18.06 x 10 22 /cm 3 ). After substitution one obtains 
hw c /2k B = 1.3K what is close to the critical temperature T c = 1.2K for Al. It allows to formulate 
the following hypothesis concerning the parametrization of the BCS model: 

hu c ~ g ~ A(0) ~ 2k B T c . (10) 

Notice that fllOj) agrees with the suggestion hw c ~ k B T c in the original BCS paper, as well as with 
the BCS relations ([J]). On the other hand the assumption u c = u>d (Debye frequency) usually 
made in the modern literature implies 

hu c = hu D » g » A(0) ~ 2k B T c . (11) 

The choice Hlu c ~ g allows to replace the electron kinetic energy in the BCS Hamiltonian by a 
constant Ep what leads to a much simpler exactly solvable model which predicts the relation 
9 = A(0) = 2k B T c [sj. 

Another consequence of ( JTUj) is the following plausible picture of electrons recombination into 
Cooper pairs. For the temperatures T > T c only about 4k B TN(0) ll| electrons are thermally 
excited and contribute to specific heat and conductivity; the rest is "frozen in the Dirac see". At 
the temperature T c those 4k B T c N(0) = 2fru) c N(0) = K electrons begin to recombine such that at 
zero temperature all of them form K/2 Cooper pairs while the rest of electrons remain frozen all 
the time. 

The next consequence of (fTUj) is presented on FIG. 2 which shows the relations between the 
coupling constant g and the critical temperature T c under the hypothesis (TIDl and (fill) , respectively. 




FIG. 2. Critical temperature as a function of the coupling constant (in Debye energy units Ed = Hlod) 
for ksT c = 1.13 Ed /2sinh(i?£) / 1 g) (solid line) and ksT c = g/2 (dashed line). 

Obviously, (TTUj) allows to obtain quite high critical temperatures at the cost of moderate increase 
of the electron-phonon coupling, while (TTTj) strongly suppresses the increase of T c . For example, to 
obtain the critical temperatures 1.2K(Al) and 9.8K(Nb) one needs, assuming ( ITT]) , the coupling 
constant g/ks equal to 72K and 80fT, respectively. The same values of g, assuming f lTOj) . yield 
critical temperatures 36K and AOK. Therefore, (fTOl) suggests that the phonon-mediated pairing 
mechanism could be sufficient to explain high-temperature superconductivity, at least for a certain 
class of such materials [12| . 

In conclusion, it has been shown that the experimental results for a small Josephson junc- 
tion loop interpreted in the light of the presented model do not only support the original BCS 
parametrization, what can be relevant for the theory of high-temperature superconductivity, but 
also change the physical picture of a flux qubit itself. Instead of the standard interpretation in 
terms of two macroscopic quantum states corresponding to the motion of millions of Cooper pairs, 
here the qubit's basis consists of the ground state and the state which differs from the ground one 
by a single excited pair, only. Obviously, more experimental evidence concerning flux qubits of 
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different sizes and made of different materials is needed as well as the first principle theoretical 
justification of the hypothesis (TTUj) . 
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